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ABSTRACT 

The early reionization of the intergalactic medium, which is favored from the 
WMAP temperature-polarization cross-correlations, contests the validity of the stan¬ 
dard scenario of structure formation in the cold dark matter cosmogony. It is difficult 
to achieve early enough star formation without rather extreme assumptions such as 
very high escape fraction of ionizing photons from proto-galaxies or a top-heavy initial 
mass function. Here we propose an alternative scenario that is additional fluctuations 
on small scales induced by primordial magnetic fields trigger the early structure forma¬ 
tion. We found that ionizing photons from Population III stars formed in dark haloes 
can easily reionize the universe by 2 ~ 15 if the strength of primordial magnetic fields 
is between 0.7-1.5 x 10 -9 Gauss. 

Key words: stars: formation - galaxies: formation - large-scale structure of universe: 
magnetic fields - cosmology: theory 


1 INTRODUCTION 


The reionization process of the intergalactic medium (IGM) is one of the major remaining problems in modern cosmology. 
From the Gu nn-Pe terson test of QSO absorption lines, it is known that the vast majority of IGM is ionized by z ~ 6 
feecker et aljEfflOll : iFan et aDbopj ) . The recent measurement of the cosmic microwave background (CM B) temperature and 
polarization cross correlations by WMAP implies that the optical depth of the universe is about 0.17 JSnergel et alJ 1200.4 
iKotmt, et aljEfifi^ ). This result favors the early reionization scenario: the reionization process occurs at z = 15 - 20 and the 
reionization s ources are first stars, unlike quasars and galaxies which are known as the reionization sources previously (for 
the details see lboeb <G Barka.nalkoOll l. 

The early reionization process by the stellar sources has been studied in detail after WMAP llGenl2f)04lFukugita fc Kawasaki 
1200.4 Idardi et alJl20o4 ISomerville fc Livioll20f)4 lHaiman fc Holdeill2004) . In these works, cold dark matter cosmogony with 
WMAP parameters is employed. What they found was it is difficult to get r = 0.17 if the standard Salpeter initial mass 
function (IMF) is adapted. To have early enough reionization, one needs to assume almost 100% escape fraction of ionizing 
photons from proto-galaxies, or introduce a top-heavy IMF. Heavy stars may form in the early universe induced by the H 2 
molecular cooling while it is still little known about the IMF in the early universe. 

An alternative scenario to realize early reionization is to enhance the amplitude of the CDM power spectrum on very 
small scales. Such enhancement makes the dark haloes form earlier. Accordingly the star formation process starts early 
enough t o achieve r = 0.17 . Th e observa tional data from redshift surveys of galaxies such as 2dF and SDSS, and Ly- 
Q clouds llSnergel et a,l. 1 1200.4 ISelia.k et al H) strongly constrain the amplitude of the density fluctuations in the scales 
larger than 1 Mpc. However, the amplitude in the scales which are relevant to the first star formation, 0.01 - 0.1 Mpc, is 
still unclear. Therefore there is still room for considering additional power in the power spectrum on very small scales. For 
example, models with initially running spectral index n > 1 on small scales, additional fluctuations from t he isocurvature 
modes, or n on- Gaussi an statistics are considered in the context of early reionization llGhen et, alJl20o4 lAvelino fc Liddlel2004l : 
ISugivama et alJl2004ll . 

If there exist strong enough primordial magnetic fields at the reionization epoch, these magnetic fields produce additional 
density fluctuations of baryons by the Lorentz force. The magnetic tension is more effective on small scales where the entan¬ 
glements of magnetic fields are larger. Therefore we expect to have additional power in the density power spectr um o n small 
scales as is the case of isocurvature perturbations. Structure formation by magnetic fields was first discussed bv IWassermanl 
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(1978). More detailed a nalysis was car ried by Kim_et_aD jl99|j) and the infl uence on the formation of large scale structure was 
recently estimated bv lOonal fc SethU ll200.1l) . Tset,hi fc Subramanianl ||200 ItI 1 pointed out that nanoGauss magnetic fields can 
induce early structure formation and may have the potential to achieve the early reionization implied by the WMAP team. 
They also studied reionization of IGM induced by the dissipation of magnetic energy. In this paper, we thoroughly investigate 
the role of primordial magnetic fields on the early reionization process. We concentrate on the early structure formation due 
to the additional power spectrum generated by magnetic fields. 

It is known that there exist magnetic fields of several /rGauss in most of the galaxies and the clusters of galaxies while 
the origin of these magnetic fields is still uncertain. Coherence lengths of these magnetic fields are typically 100 kpc - 1 Mpc 
jKronber ellll994ll . 

Perhaps small seeds of the magnetic fields are produced inside astronomical objects such as stars and AGNs due to the 
Bierman n battery me chanism. Although the resultant magnetic fields ar e very weak, those may be amplified by the dynamo 
process iZeldovich et al.lil98.il : for a comprehensive review see IWidrowf 12002 1. Eventually these magnetic fields are spread 
by Super nova winds or AGN jets into IGM. However , for achieving observed values in cluster s of galaxies and high red shift 
galax ies iKronberg et a j] ll992 : iKim et, aklll99fll . Il99lf) . there are difficulties in dynamo theory ferandenburg fc Subramanianl 


An alternative to the dynamo scenario is the generation of magnetic fields in the early universe, which we consider in this 
paper. Magnetic fields can be formed either due to the bubble collisions during the cosmic phase transition such as QCD or 
electroweak phase tr ansitions, or due t o the break of the conformal invariance in the Maxwell theory during the inflation. For 
a detailed review, see iGiovanninill2004l These magnetic fields are formed early enough to make influence on the first structure 
formation in the universe and the reionization process. 


The primordial magnetic fields are constrained bv Bit 

t-Bang nucleosynthesis and CMB 1 

Cheneet al. 1996:; Kcrnan et al. 

199ft .Tedamzik et al. 2000: iMack et al. 2002!; Lewis 20041: 

Yamazaki et al. 2005; iTashiro et al. 

200JJ). The upper limit of the 


comoving amplitude of magnetic fields is ~ 10 nGauss at the 1 Mpc scale. 

Throughout this paper we use the cosmological parameters measured by the WMAP teams,: the Hubble co nstant (in the 
unit of 100 kms _1 Mpc _1 ) h = 0.71, the matter density ratio = 0.27 and the baryon density ratio fib = 0.044 JSoergel et alJ 

l200.il . 


2 THE DENSITY FLUCTUATIONS GENERATED BY THE PRIMORDIAL MAGNETIC FIELDS 

Primordial magnetic fields produce the density perturb ations after recombination ( Wasserma nil 978trKim et, a, ill 99fil : ISuhra,ma,nia,n Rr. Barrow! 
1998 : i Oopal fc Sethil2Q*oii : ISethi fc SubramaniarfePOBlf by the Lorentz force. The generated perturbations grow gravitationally 
and end up collapsing to form structure. The density fluctuation evolution with the primordial magnetic fields is governed by 
following equations, 


d 2 S h 
dt 2 


— 2— — 7 ^ - b 47rG(pb<5b T" Pdm5dm) + S (t , x) , 


S(t,x) 


V- [(V x Bo(x)) x Bo(x)] 
47rpboa 3 (f) 


(1) 

( 2 ) 


<9 2 <5dm 
dt 2 


2 a d6dm 
a dt 


T 47rG(pb^b T Pdm^dm), 


( 3 ) 


where a(t), 5 and p are the scale factor, the density perturbations and the energy density, and the subscripts b, dm and 0 
denote the baryon and dark matter components, and the present (comoving) value. The dot represents the time derivative. 

In Eq. Q. we combine equations for ionized and neutral baryons. In the early universe, the interaction between ions and 
neutrals is so efficient that they are tightly coupled and move together. The interaction rate is proportional to the relative 
velocity between ions and neutrals. Since we can assume the thermal velocity of baryons as the relative velocity if there is 
no magnetic fields, the interaction rate, which is also proportional to the residual ionization fraction, becomes shorter as the 
universe expands. Eventually, the time scale of the interaction (inverse of the interaction rate) becomes longer than the Hubble 
time at z « 10 for the value of the ionization fraction x e « 10 -4 . Since then, we need to treat ions and neutrals separately. 
However, if there exists magnetic fields, this is not the case because the velocity of ions (and the relative velocity) becomes 
almost the Alfven velocity va = B /y / 47rpb*7 which is 1000 times larger than the thermal velocity at z ~ 10. Accordingly the 
time scale of the interaction becomes 1000 times shorter and never exceeds the Hubble time before the reionization. 

In order to solve these equations, we define total matter density and its perturbation p m and <5 m as 


Pm — (Pb T Pdm); 
x _ [pbSb + Pd m 5dm) 

Om = -• 


( 4 ) 


( 5 ) 
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We can write the equations for <5b and <$ m from Eqs. 0 and 0 as 

d dfib | A S -1 r I q(+ \ 

-oTo- = -2-—- +47T Gpmfim + S(t,x), 

at 2 a ot 

d ^ ^ , Pb cU 

r^.,^ — 2 + 47rGp m ^m + o(t, xj. 

at 2 a at Pm 

The solution of Eq. q can be acquired by the Green’s function method, 

i-t 


<5 m = A(x)Di(t) + B(x)D 2 (t) - £>i(t) 


11b 
- -11 


df‘ 




dt 


' i j 


W(t') 


w(t') 1 ‘" vv n n 

where -Di(f) and D 2 (t) are the homogeneous solutions of Eq. 0 , W is the Wronskian and is expressed as 

w(t) = Di(t)i) 2 (t) - D 2 (t)i>i(t), 


( 6 ) 

(7) 

( 8 ) 

(9) 


and ti denotes the initial time. 

The first and second terms of Eq. 0 are corresponding to the growing and the decaying mode solutions of primordial 
perturbations produced by inflation and the third and fourth terms are the ones generated by the primordial magnetic fields. 
Hereafter we only consider the growing mode solution for the former terms and describe it as 5 m p. The later two terms are 
written as <5 m M- 

Since we only treat the evolution of perturbations in the matter dominated era, D\{t) oc t 2 ^ 3 and D 2 (t) oc f -1 while they 
should be modified once th e cosmological constant te rm starts to dominate in expansion of the universe, z<0.5. Accordingly 
we obtain <5 m M as Iwassermanlll978t iKim et alJll996 1 


x fib 

OmM = — 


— (L \ 2/3 -fl\ 
To \TJ + 5VtJ 


f 2 S(fi,x). 


( 10 ) 


The density fluctuations generated by the primordial magnetic fields have the same growth rate as the primordial density 
fluctuations. 

Next, let us calculate the power spectrum of the matter density perturbations. For simplicity we assume that there is no 
correlation between the magnetic fields and the primordial perturbations. With this assumption, the matter power spectrum 
can be described as 


Pm(fc) = Bnp(fc) + P mM (fc) = (|<5 mP (fc)| 2 ) + <|5 mM (fc)| 2 >, 


( 11 ) 


where <5 m p(fc) and <5 m M(fc) are Fourier components and ( ) d enotes the ensemb le average. 

We obtain P m p(fc) by using the CMBfast code iselia.k fo Za.lda.rriagallT99 61 while PmM(fc) is calculated from Eq. m 


PmM (&) 


where 


/ n b \ 2 

( 

v 

r 9 

ft \ 2/3 ,3 ft\ 3] 

(fim) 

^4npboa 3 (ti) 

) 

10 

U) + 5 Vii/ 2 


I 2 (k), 


T(k) = (|V ■ (V x Bo(x)) x Bo(x)| 2 ). 


( 12 ) 


(13) 


To calculate I 2 (k), we need to know the power spectrum of the magnetic fields. The amplitude of the power spectrum 


Bg(k) is defined as 


(B oi (ki)Bo*,-(k2)> s 


(2^ 


<5(ki - k 2 ) Sij - 


k\ik 2 i 


B 2 0 (k). 


The nonlinear convolution Eq. m leads to 


I 2 (k) = / dki / dfiL 


Pg(fci)P 0 2 (|k-ki|) 


[2k 5 klfi + k k\ (1 — 5/i )-\-2k 3 ki/i ]. 


k — ki | 2 

For the power spectrum B 2 (k), we adopt the power law shape with sharp cutoff at k c as 
k ' n+3 


( h \ 

j-J 


Bn 


for (k < k c ), 


(14) 


(15) 


(16) 


where n is the spectral index. Here Bo is the rms amplitude of the magnetic fields in real space as (B 2 (x)} = B 2 , which 
provides the characteristic magnetic amplitude at the cutoff scale. 

The cutoff scale depends on the dissipation mecha nism of the magnetic field energy ll.Teda.mzik et a.l 1 1 998tlSubra l ma,nian fc Barrowl 
ll99f4lBaneriee fc ■TedamziklteooitlTa.shiro et ahlfeoPfilf . Before recombination, dissipation is caused by the interaction between 
baryons and photons. After recombination, the magnetic field energy in small scales is mostly dissipated by the nonlinear 
effects, i.e., direct cascade. The dissipation time at the scale l due to the direct cascade is the eddy turn-over time, l/v, where 
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v is the baryon fluid velocity. Since the interaction between baryons and photons which plays a role as a viscosity is no longer 
effective, the velocity of the baryon fluids increases and the equipartition between the magnetic fields and the kinematic energy 
is established. Hence, we assume that the fluid velocity achieves the Alfven velocity. Under this assumption, the cutoff scale 
is determined as the scale where the eddy turn-over time is equal to the Hubble time P" 1 , 

H . _i i B 0 y 1 ( n h \ 1/2 / h 


k, 


2?r- 


va 


102 


Mpc " 1 ( 


1 nGauss/ 


fJlL. V (—Y 

V 0.044/ V0.71/ ' 


(17) 


Note that this comoving cutoff scale is constant in the matter dominated epoch. In this paper, we consider the power spectrum 
with sharp cutoff below l/k c for simplicity as is shown in Eq. 11 (ill to calculate I 2 {k). With this assumption, k\ integration of 
I 2 (k) is determined by the value of integrand at k = k\. However, in reality, the direct cascade process of the magnetic fields 
from large scales to small scales likely produces a power-law tail below the cutoff scale in the power spectrum. Although most 
of the contribution for I 2 (k) still comes from the peak of the spectrum, i.e,. ki = k c as well as the sharp cutoff case, some 
corrections may be needed. 

Finally we introduce one more important scale for the evolution of density perturbations, i.e., magnetic Jeans length. 
Below this scale, the magnetic pressure gradients, which we do not t ake into account in Eq. 0 , counteract the gravitational 
growth and prevent further growth of density perturbations. iKim et akllll 99fj evaluated the magnetic Jeans scale as 


A-mj = 5—-- = 32 

Bn 


Mpc " 1 ( 


Bo rVllmW i 
1 nGauss ) V 0.27 A 0.71 


(18) 


3 REIONIZATION 

The Population III stars are cons i dered as t h e main s ources of reionization. We can estimate the number of photons from 
Populations III stars following lSomerville fc LiviofcOQ.jl In their analysis, ionizing photons are derived from the star-formation- 
rate density, which is assumed to be proportional to the time derivative of the fraction of the total mass in collapsed haloes 
F. To obtain F, the Press-Schechter prescription is employed, in which additional small-scale power in the matter power 
spectrum induced by magnetic fields plays an important role. 


3.1 The mass dispersion 


From the matter power spectrum P m (fc, f), the mass dispersion within a radius R can be written as 


cr 2 (P, t) = 4n 


dkk 2 P ul (k, t) exp (—fc 2 P 2 ) . 


(19) 


Here we adopt the Gaussian window function. 

The total mass M within R in the matter dominated epoch can be described as 


R = 18 kpc 


f M N 

1/3 ^ 

\ /On\ 

" 1/3 / h V 

l 1O 6 M 0 

/ \0.27/ 

\0.71/ 


-2/3 


( 20 ) 


We plot the evolution of the mass dispersion which is calculated from the matter power spectrum induced by magnetic 
fields in Fig. Q On large mass scales, the mass dispersion is proportional to M “ 7 ^ 3 for the case of the spectral index n— 1. 
The dependence on the spectral index is discussed below. On the other hand, on the mass scales lower than the magnetic 
Jeans mass Mmj — 1.2 x 10 9 (Bo/l nGauss) 3 M 0 , the slope of the mass dispersion becomes milder due to the sharp cutoff 
of the matter power spectrum below the magnetic Jeans scale in Fourier space. For the comparison, we also plot the mass 
dispersion without taking into account the effect of the magnetic Jeans oscillations as dotted lines. 

Since the mass scales relevant to the Population III star formation are around 10 6 - 10' M@ as is shown in the next 
subsection, we hereafter focus on these mass scales. In Fig. [5] the evolution of the mass dispersion is seen for the model 
with Bo = 0.7 nGauss and n = 1. Since the growth rates of both P m p and P m M are f 4 / 3 oc 1/(1 + z ) 2 , the mass dispersion 
a 2 oc 1/(1 + z ) 2 too. In this figure, we plot the contributions from P m p and P m M separately at z = 20. It is shown that below 
M = 5.0 x 10 6 Mq, the contribution from magnetic fields, i.e., P m M is dominated. This mass scale is independent on redshift. 
Accordingly we expect to have early structure formation induced by P m M- 

We show the mass dispersion for models with different magnetic field amplitudes at z = 20 in Fig. [3] The contributions 
from magnetic fields are plotted as the solid lines. The dotted line represents the contribution from the usual primordial 
perturbations. We can analytically estimate Bo dependence from Eq. m as a 2 oc I 2 oc Bq which is consistent with numerical 
results. It is shown that the effect of magnetic fields on early structure formation or reionization is significant if Bo > 0.6nGauss 
since there appears extra-power above the smallest collapsed haloes to form Population III stars, i.e., 10 6 - 10 7 Mq. As the 
magnetic field strength becomes stronger, the magnetic Jeans scale shifts to the large scale. Accordingly, the slope of the mass 
dispersion on the relevant mass scales becomes flatter as is shown in Fig. [3 
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Figure 1. The redshift evolution of the mass dispersion of the density perturbation generated by the magnetic fields with Bo = 0.7 
nGauss and n = 1. The solid lines are the mass dispersion at z - 15, z = 20 and z = 40 (from top to bottom). The dotted lines are the 
mass dispersion without considering the effect of the magnetic Jeans oscillations. 



Figure 2. The redshift evolution of the mass dispersion for the model with Bq = 0.7 nGauss and n = 1. The solid lines represent the 
total mass dispersion at z = 15, z = 20 and z = 40 (from top to bottom). The dotted and dashed lines are the contributions from P m p 
and P m m, respectively, at z = 20. 


The mass dispersions for various spectral indices at z = 20 are plotted in Fig. HI The dependence on n can be e xplained as 
follows. In the limit of k/k c 1, I 2 (k) ~ ak 2n+7 + /3/c 4 where a and f3 are coefficients which depend on n and Bo (iKim et alJ 
Il99d iGooal fc Sethill2003l *l. The former term dominates if n < —1.5, while the later one dominates for n > —1.5. From the 
definition Eq. 111 ) 1 . the mass dispersion er 2 oc fc 3 P m (fc) oc k 3 I 2 (k). Accordingly a 2 oc oc M~ 7 ^ 3 for n > —1.5 and 

^. 3 +( 2 ra+ 7 ) x JV/-( 2 n +io)/3 f or n < — 1 . 5 . it i s interesting that the slope of cr 2 does not depend on n for n > —1.5. Note that 
the slope becomes milder in the small mass around the Jeans scale regardless of the value of n as mentioned above. 



M [M 0 ] 


Figure 3. Mass dispersions of models with different strengths of the magnetic fields. The solid lines show the contributions from the 
magnetic fields of 1.0 nGauss, 0.8 nGauss, 0.7 nGauss and 0.6 nGauss (from top to bottom) at 2 = 20. 
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M [M 0 ] 


Figure 4. Mass dispersions of models with different spectral indices for the power spectra of the magnetic fields at z = 20. We fix the 
strength of the magnetic fields as 0.7 nGauss. The dotted, dotted-dashed, solid, dashed lines (from top to bottom) are the models with 
n = —2, —1, 1 and 2, respectively. The thin dotted line is the contribution from the primordial perturbations. 


3.2 The ionization photon number 


We assume that the Population III stars are formed in collapsed haloes with masses larger than M cr it and lower than M v n. 
Here we adopt M cr it = 1.0 x 10 6 and Af v i r = M(T V u = 10 4 K) where Af(T v j r ) is the virial mass with temperature 

T v j r . The virial mass M(T V1V ) evolves as (1 + z) _1 for given temperature T v i r . At 2 = 15, the virial mass with T v u = 10 4 is 
M (T vir = 10 4 ) = 5.2 x 10 6 A4>. 

Following the Press-Schechter prescription, we can derive the fraction of the collapse haloes with masses larger than M 
at time t as 


F(> M,t) = 1- 


dx- 


V2 


t ) 


exp 


2 a 2 (M,t) 


( 21 ) 


where cr 2 (M,t) is obtained by substituting Eq. 12()t into Eq. IPl . 

From F, the global star-formation-rate density of the Population III stars can be calculated as 

p*(f) = e *Pb^- [F(> M cri t, t) - F(> M vir ,f)], 


( 22 ) 


where e* is the efficiency of conversion of gas into stars and we adopt e* = 0.001 lIYoshida et a,l.1120031 1. 

Now we assume that the lifetime of Population III stars is A t = 3 x 1 0 s yr. and the production rate of ionizing photons by 
Population III stars is IV 7i o = 1-6 x 10 48 photons s~ 1 M^ 1 (Bromm et al ] 2ljoij). Accordingly we obtain the total production 
rate of ionizing photons as 


n 7 (t) = dt' iV 7> o0(Af — (f — f , ))p*(f , ) 

= TV 7 ,o [p* (f) - p, (t - At)], 


(23) 


where 6 {t) is a step function which is 6 (t) = 1 for t > 0 and 9{t ) = 0 for t < 0. 

The cumulative number of photons per H atom is a good indicator of the ionization degree of IGM. It is kno wn th a t about 
20 cumulative p hoto ns per H atom are requited in order to achieve a volume weighted ionization of 99 percent llHaiman et al] 
2001; Sok asia.n et a 1 .1120021 . 1200417 From Eq. 1231 . the cumulative photons per H atom at time t is represented as 


^cumul 


n h 


(*) 


P»tp 

Pb 



dt'h-y(t') 


= 1V 7 ,[ dt 1 lp*(t') — p*(t' — At)l , (24) 

Pb Jo 

where hh is the Hydrogen number density, p is the mean molecular weight, and m p is the proton mass. 

We show cumulative photons from the Population III stars per H atom as a function of redshift. In Fig.^1 the cumulative 
photons are plotted for different values of the magnetic strength Bo. Here we fix n = 1. It is found that more photons are 
produced as the magnetic field strength becomes stronger since the amplitude of the density perturbations depends on the 
strenth. We can conclude that the universe is reionized early enough to be consistent with WMAP data if Bo is larger than 
0.7 nGauss. Remember that 20 cumulative photons are needed to achieve a 99 percent volume weighted ionization. 

If the magnetic field strength is larger than InGauss, however, it is shown that the cumulative photons become smaller 
as we increase the strength since the magnetic Jeans scale shifts to the larger scale which makes Population III stars difficult 
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Figure 5. The cumulative photons from the Population III stars per H atom for the models with the magnetic fields. The spectral 
index of the magnetic fields is fixed as n = 1. The strengths of the magnetic fields are 2.0, 1.5, 0.8, 0.7 and 0.6 nGauss for the dashed- 
dotted-dotted, dashed-dotted, dotted, solid and dashed lines, respectively. The thin dotted line is the case with no magnetic fields. 



Figure 6. The cumulative photons from the Population III stars per H atom for the models with the fixed strength of magnetic fields 
0.7 nGauss and various spectral indices. The dotted, dotted-dashed, solid, dashed lines (from top to bottom) are the models with n = —2, 
— 1, 1 and 2, respectively. It is shown that dependence on n of the cumulative photons is very weak. The thin dotted line is for the case 
of no magnetic fields. 


to be formed. We found that the magnetic fields larger than 1.5 nGauss cannot generate large enough number of cumulative 
photons to induce early reionization. 

Fig. El shows the dependence of the cumulative photons on n. It is found that there is little dependence if n > —1.5. 
Therefore the requirement Bo > 0.7 nGauss for early reionization is robust regardless of the power law index of the magnetic 
fields. 


4 CONCLUSION 

In this paper we investigate the role of the additional density perturbations generated by the primordial magnetic fields on 
the reionization process in the early universe. These additional density perturbations may trigger the early structure and star 
formation. Employing a simple analytic recipe, we estimate the number of ionizing photons emitted from the Population III 
stars. We found that the reionization process almost completes by z ~ 15 if the strength of primordial magnetic fields is 
larger than 0.7 nGauss and less than 1.5 nGauss. Note that we adopt the Gaussian window function to calculate the mass 
dispersion (see Eq. ITTil'i. Different choice of the window function alters the mass dispersion at the magnetic Jeans scale. 
Accordingly the magnetic field strength to be required for the early enough reionization is also changed. If we employ the 
sharp-fc window function, for example, the strength should be between 0.5 and 0.7 nGauss. 

Such magnetic fields are not yet ruled out from current observations, i.e., BBN, CMB temperature anisotropies and 
polarization, and Faraday rotation of polarized lights from radio sources. Although the formation process of such primor- 
dia l magnetic fie ld s is still uncertain, magnetic fields may be naturally generated during the cosmological phase transition 
feaneriee fc JedamziklEoCMlh 

The extra-power of the matter power spectrum will be directly probed by future observations such as the fluctuations of 
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the Hydrogen 21cm line llLoeb & Zaldarriagal2Qo3l and the substructure of lensing haloes (iDalal fc Kochanekl2002() . Moreover, 
the thermal diffusion process of primordial magnetic fields may cause ionization of IGM even before Popu lation III stars. A 
measurem ent of fluctuations of the 21cm line will be a powerful tool to investigate such pre-reionization (iTashiro fc Sueivamal 
lin preparation!). 
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